Muon Spin Relaxation Studies of Superconductivity 
in a Crystalline Array of Weakly Coupled Metal Nanoparticles 
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We report Muon Spin Relaxation studies in weak transverse fields of the superconductivity in 
the metal cluster compound, Ga84[N(SiMe3)2]20-Li6Br2(thf)20-2toluene. The temperature and field 
dependence of the muon spin relaxation rate and Knight shift clearly evidence type II bulk super- 
conductivity below T c ~ 7.8 K, with B c i w 0.06 T, B c2 ~ 0.26 T, k ~ 2 and weak flux pinning. The 
data are well described by the s-wave BCS model with weak electron-phonon coupling in the clean 
limit. A qualitative explanation for the conduction mechanism in this novel type of narrow band 
superconductor is presented. 

PACS numbers: 76.75.+i, 74.78.Na, 74.10.+V 



The chemical synthesis of molecular metal cluster com- 
pounds presents an attractive bottom-up route for the 
generation of self-organized nanostructures composed of 
3D ordered arrays of identical metal nanoparticles em- 
bedded in a dielectric matrix. Until recently, such cluster 
solids were always found electrically insulating. On the 
other hand, the strong similarity with (super) conducting 
molecular crystals, as the alkali-doped fullerenes (Ceo), 
suggests that in principle metal cluster compounds could 
also display metallic conductivity (and even supercon- 
ductivity) due to intcrmolecular charge transfer. We 
have indeed recently obtained compelling evidence from 
69,71 Ga-NMR [l| and magnetization measurements for 
the occurrence of band-type conductivity in crystalline 
ordered Gas4 cluster compounds, composed of arrays of 
giant Gag4 cluster molecules that display mixed-valence 
properties. In addition, bulk type II superconductivity 
was observed below a transition temperature T c w 7.5 K, 
much higher in fact than known for bulk a-Ga metal 
(T c « 1 K). This material may thus represent a first ex- 
perimental realization of a theoretical model advanced by 
Friedel in 1992 |2j, who predicted that for a crystalline 
array of identical metal nanoparticles, a very weak inter- 
particle charge transfer can still yield superconductivity 
with a relatively high T c value. 

From the NMR study, a number of results were ob- 
tained. First, the field and temperature dependence of 
the Knight shift can be well fitted to BCS theory for 
weak electron-phonon coupling, a rather surprising re- 
sult since strong correlations would be expected for such 
narrow band conductors as the present. Second, a very 
strong sample dependence of the second critical field B c2 
was revealed. As compared to B c2 « 13.8 T reported 
previously [jj , we found that for most of the samples the 
minimum field of about 2 T required for the NMR sig- 
nal (due to signal-to-noise ratio) was already sufficient to 
suppress the superconductivity completely. Only in a few 



samples superconductivity could be detected by NMR up 
to about 5 T. By contrast, magnetization studies showed 
diamagnetic signals for all samples below similar values 
of T c —7-8 K, but with apparent B c2 values as low as 
0.3 T. Although an explanation of this variation (by a 
factor 50!) in terms of lattice defects — associated with 
orientational disorder in the molecular crystal lattice — 
was proposed, definite proof for the occurrence of bulk 
superconductivity in the low B c2 samples seemed highly 
desirable. In order to check this assumption, muon-Spin- 
Relaxation (/iSR) experiments in small transverse fields 
(Btf) have therefore been performed. As discussed be- 
low, they provide unambiguous proof for bulk type II su- 
perconductivity, with a Ginzburg-Landau (GL) param- 
eter k ~ 2, and are likewise in agreement with weak- 
coupling BCS model predictions. A qualitative discus- 
sion of the nature of the conduction mechanism in terms 
of intcrmolecular charge fluctuations is presented at the 
end of the paper. 

The crystal structure of the giant molecular cluster 
compounds Gas4, which is short for Gas4[N(SiMe3)2]20- 
-LieBr2(thf)2o-2toluene, is fully described in Q|. The 
/iSR experiments were performed on the sample labelled 
S3 in the previous work Q], which showed an excellent 
crystallinity in x-ray experiments with typical dimen- 
sions of the crystallites of 0.1 mm. Magnetization mea- 
surements in a commercial SQUID magnetometer evi- 
dence a SC transition at T c = 7.75(5) K in zero external 
field. A T — extrapolation yields an upper critical 
field B c2 w 0.26 T. Using the expression B c2 ps $ /27t£ 2 , 
where $o is the flux quantum, gives a coherence length 
£ ~ 35 nm much larger than the inter-cluster distance, of 
order 2 nm. The extrapolated thermodynamics! critical 
field at T — was found as B c \ w 60 mT. Using the gen- 
eral expression given by Brandt for an ideal triangular 
GL vortex lattice, relating B c2 , B c \ and the GL parame- 
ter k p| , we find n ~ 2 and a London penetration depth 
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A ~ 70 nm. This sample, which is consequently close 
to the "clean" SC limit, indeed displays the largest con- 
ducting fraction reached so far in a Ga 8 4 sample (90%), 
as found in the 69,71 Ga NMR experiments performed on 
several batches Q . 

Transverse field /iSR is a powerful technique to probe 
type II superconductivity locally and can be per- 
formed in any arbitrary field close to zero, contrary to 
NMR. By implanting 100% spin polarized muons in 
the material, the local field they experience can be mea- 
sured through their decay into, among others, a positron, 
in a typical time window of 0.1-10 /is. We performed our 
experiments on the 7rM3 beam line on the GPS spec- 
trometer in PSI, Villigen, Switzerland, in the temper- 
ature range 2 K < T < 10K, using magnetic fields 
< Btf < 0.4 T. The Ga 8 4 sample (30 mg) is extremely 
air sensitive and must be kept in a toluene solution with a 
sample/toluene mass ratio of order 80/20(±10) to avoid 
loss of crystal solvent. We therefore sealed it in an Al 
sample holder [2^, using an Al thickness of 0.35 mm 
between the incident /i + beam and the sample as a mod- 
erator for the muons. From the geometry our sample, 
a pellet of thickness 0.55 mm and diameter 7 mm, and 
the calibrated beam distribution, the n + fraction stopped 
in the Al sample holder (plus the non conducting frac- 
tion of Ga 8 4) called "background" in the following, is of 
order 25%. This will appear as a temperature- and field- 
independent signal in the whole temperature range. The 
/i + fraction stopped in toluene is of order 17% and does 
not contribute to the measured signal, as seen experi- 
mentally by a reduction of the expected intensity by this 
amount. We attribute this reduction to the formation 
of muonium in this solvent. Given the large fraction of 
"lost" muons (~ 42%) we took data with good statistics, 
with typically 40 x 10 6 [i + per point. 

Figure ^ shows the time dependence of the fi + polar- 
ization P(t) along the axis defined by Btf, for two tem- 
peratures and -Btf = 60 mT > B cl , in a Field Cooled 
(FC) experiment. A clear difference can be seen between 
the data measured at T = 10 K > T c and T = 2 K < T c . 
Above T c , the local field (B M ) distribution is uniform in 
the sample and centered at the corresponding Larmor 
frequency jB^/2-rr, where 7 is the /z + gyromagnetic ra- 
tio. In real time, the relaxation of P(t) is therefore very 
slow. Below T c , a flux-line lattice (FLL) is created in the 
mixed state, resulting in a larger distribution of B^, i.e, 
a faster relaxation of P(t). 

In order to get a first qualitative impression of the 
behavior of the various relevant parameters, a single 
damped cosine was fitted to the raw data. This model- 
independent analysis already clearly shows a transition 
at T c w 6.4 and 5.9 K, for Btf — 60 and 90 mT (cor- 
responding to B C 2 found in magnetization measurements 
at these temperatures), in a FC experiment, in both the 
relaxation rate (not shown) and the precession frequency, 
Hot [Fig^inset)]. In a second step, the background con- 
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FIG. 1: Polarization P(t) of the /i + spins, in a FC experi- 
ment (Htf = 60 mT). For clarity, only one of the two parts 
(imaginary or real) is presented for each T, in a rotating frame 
at a frequency 7 x 55 mT. The continuous lines are fits with 
Eq.0 The dotted and dashed lines represent the envelope of 
the "background" and the total signal, respectively. 

tribution has to be evaluated to get more quantitative 
information. To fit the /i + polarization, we used the phe- 
nomenological function 

P(t) = xexp(-[cr Ga i] Q - er 2 d i 2 )cos(w Ga i + </>) 

+ (1 - x) exp(-crg gd £ 2 ) cos(u hgd t + 0) . (1) 

The first term accounts for the Ga 8 4 signal. We 
point out that these fj, + are likely stopped in the 
Ga 8 4[N(SiMe3)2]2cT clusters, the other locations in the 
crystal corresponding to regions having a positive charge. 
The second term represents the background. The relax- 
ation rates cr n( j and <7b g d are due to the nuclear dipolar 
fields on the respective sites . We first determined the 
T- and -Btf -independent parameters by comparing all 
the data. The fits yield a fraction x « 0.7 of /i + stopping 
in Ga 8 4, cr n d ~ 0.06 /xs _1 and <7b g d ~ 0.04 M s_ ii which 
are typical values for nuclear relaxation rates 0, S @- 
Although a rough approximation, a Gaussian (a = 2) is 
usually suited to fit the FLL field distribution in pow- 
der samples 0, . This implies that the second moment 
of the field distribution, (AB 2 ) 1 / 2 , related to the pen- 
etration depth A, is directly proportional to the fi + re- 
laxation rate due to the FLL distribution (<7Ga)- In our 
case a Gaussian cannot account for the decay of the first 
P(t) oscillations at low-T (Fig. ^) and we find instead a 
T- and ifrF-hidependent a ~ 0.5. This suggests a more 
complex FLL distribution |lp| in this cluster compound. 
A complication of this a value is that the correspond- 
ing field distribution does not converge unless a cutoff 
(which could not be measured within the experimental 
resolution) is introduced in the spectrum. To proceed, 
we therefore assume the proportionality <7Q a oc (AB 2 ) to 
be still obeyed. We are then finally left with only two 





FIG. 2: T-dependence of Kg& — Kbga, for two values of Btf. 
Lines are BCS fits with a s-wave symmetry of the SC gap. 
"W" and "S" are for "weak" and "strong" electron-phonon 
coupling, "C" and "D" for "clean" and "dirty" limits. Insets: 
precession frequency u to t of the spins obtained with a sin- 
gle damped cosine fit, for Htf = 60 mT. The line shows the 
average value for T > T c . 



free parameters, cca and vq & . 

The T-dependence of the average shift Kq & — (uJGa — 
iBtf)/iBtf, is displayed in Fig.|2|with the background 
shift i^bgd subtracted. This (muon) Knight shift is usu- 
ally the sum of a T- and £?TF-independent orbital shift 
-ftTorb and a spin part K$, proportional to the s-electron 
susceptibility [and to the electronic density of state at 
the Fermi level, D(Ep)]. In the normal state (n), we find 
a T- and £>TF-independent value (as expected in met- 
als) of i^Ga,n ~ -f^bgd ~ —0.17%. The negative sign is 
due to the hyperfine coupling of the /i + with the carri- 
ers, which is known to be strongly sample dependent [llj 
|23j. The present value is comparable to those found in 
several heavy fermion compounds fr2l . 

Below T c , the average precession frequency (inset of 
Fig. |2l is seen to increase continuously down to T — > 
for both values of Btf- Since the Knight shift is nega- 
tive, this implies a decrease of its absolute value. In the 
mixed state, a possible additional shift (-Kd) due to the 
demagnetization and Lorentz fields would reduce the lo- 
cal field at the muon site. Therefore, it cannot account 
for an increase of vtot and hence the variation of the shift 
can be attributed to a decrease of the spin part \Ks\ due 
to spin-singlet-pairing of the quasiparticles below T c [24| . 

Figure |3 shows that K S (T) is well fitted by a BCS 
model with s-wave symmetry of the SC gap and weak 
electron-phonon coupling (i.e., with a SC gap Ao « 
1.76fc B T c ) 0O]1. The decrease of the SC carrier density 
due to the increase of Btf is directly observed through 
a non-zero shift at T — > 0. It is also clear that the reduc- 
tion of | .Kg | is larger in the "low" field measurement, as 
expected from a larger fraction of Cooper pairs |l6| . The 



FIG. 3: T-dependence of /(6)<r G a oc A . The lines are BCS 
fits (WC). Inset: a Ga (T) from Eq.[T] 



strong electron-phonon coupling and the dirty limit cases 
(the latter being unlikely from the magnetization data) 
|l7| are also presented in Fig. [3 Their steeper variation 
below T c clearly does not agree with the data. 

The T-dependence of oca is presented in the inset of 
Fig. |21 The reduction of T c and an increase of the inter- 
nal field homogeneity (i.e., a decrease of (TQ a ), is clearly 
observed when Btf is increased. In the extensively stud- 
ied cuprates and heavy fermions, the conditions k 3> 1 
and -Btf •C B C 2 imply that the ratio OGa/A~ 2 is roughly 
T- and Btf -independent 0,0- Here, the artificial re- 
duction of (TQa due to the overlap of the vortices near 
B C 2 has to be taken into account. In the framework of 
the GL theory, Brandt derived the quantity \ 2 (AB 2 .) 1 / 2 
as a function of k and b, where b is the average internal 
field divided by B c2 {T) 0. Calling this quantity /(&), 
the function A~ 2 oc f(b)<jQ a is presented in Fig. |3| tak- 
ing k — 2. The qualitative variation of A~ 2 (T) is very 
similar to the T-dependence of Kg, as expected in the 
BCS model in small fields 0, 0] , and as shown by the 
fit with weak electron-phonon coupling in the clean limit 
(Fig. |3 0. 

We finally discuss the field dependence of oca at low-T. 
We measured aca from the pure diamagnetic state up to 
the normal state by increasing the field from a Zero Field 
Cooled (ZFC) condition at 2 K, from up to B T f > B c2 . 
The results of the fits with Eq. ^are presented in Fig. 0] 
For small fields (Btf ^ B c i), (TQa continuously increases 
with Btf- A strong inhomogeneity of the electron den- 
sity of state in Gag4, as well as the powder nature of 
the sample as suggested in 0, may lead to some FLL 
bending around SC grains and an inhomogeneity respon- 
sible for a non-zero CQa- When Btf ^ B ci , UQa reflects 
both the FLL organization || and the field dependence 
of A~ 2 . For high fields {B T f > B c2 - 0.2 T at 2 K), the 
superconductivity is destroyed and an homogeneous field 
(&Ga = 0) is recovered. The ratio A/A(0) oc [/(6)<TGa] -2 
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is presented in Fig. 0] as a function of b. In a BCS model 
with s-wave symmetry, the penetration depth A is ex- 
pected to be weakly field-dependent 0] , with a roughly 
linear variation of oca- However, a linear variation of A, 
quantitatively comparable to the type II superconduc- 
tor NbSe2 0|, is measured and is consistent with a d- 
wave symmetry of the SC gap (continuous lines) [fl fli^]. 
Moreover, the predicted ac a variation, with a constant A, 
does not fit the data (dashed line) . The theoretical com- 
putations are usually performed in an equilibrium FLL, 
which does not really correspond to our case. Indeed, 
the pinning of the vortices increases the local field dis- 
tribution, as seen by the usual 0, |2£| reduction of o-Qa 
in the FC case (which also reduces the sensitivity of the 
experiment). However, the shape of the variation does 
not seem to be changed. 

Finally, we address the nature of the conduction mech- 
anism in Ga84- In the Friedel model, referred to in the 
introduction, conduction would occur in a narrow band 
originating from a molecular energy level near Ep of the 
Ga84 cluster, with a width proportional to the interclus- 
ter , i.e. anion-anion charge transfer t. In view of the 
ligand shells around the Gas4 cluster cores, t is however 
expected to be quite small, of order 0.01 eV. In a Hub- 
bard model approach, t would have to compete with the 
on-site Coulomb interaction (Hubbard-t/) of the anion, 
estimated at about 0.5 eV (three times smaller than for 
Ceo)- In such a narrow-band system, strong electron- 
electron correlations are expected whereas the NMR and 
/iSR experiments both indicate nearly free electron be- 
havior (Korringa constants, Knight shifts). It appears 
relevant, therefore, to consider other possible charge fluc- 
tuations. 

In analogy with the well-known models proposed for 




B TF (kG) 

FIG. 4: HTF-dependence of UGa, in ZFC and FC conditions. 
The pinning accounts for the difference. Inset: A(6)/A(0). 
The lines are computations with k = 2, a constant A (s-wave, 
dashed lines) and A/A(0) oc 1 + 3.1(5)6 (d-wave, solid lines). 
We note that b = when B TF < B cl 0. 



the transition metal compounds 21], one such possibil- 
ity could be an electron transfer from the cluster anion 
to the Li-cation, in combination with conduction in the 
cation band derived from the cation-cation overlap. Pre- 
liminary estimates, based on the cation ionization energy, 
the anion electron affinity and the Madelung energies in- 
volved, indicate the anion-cation charge transfer process 
to be almost energetically neutral, i.e. could in fact be 
quite small, smaller perhaps than the width of the cation 
band. The latter is difficult to estimate but could be 
larger than t. In this scenario the itinerant carriers would 
reside mainly in the cation band, which could explain 
why in the measured Ga-NMR spectra the different Ga 
sites in the cluster are clearly resolved, as in an insulating 
Ga compound. Clearly, more sophisticated calculations 
are needed to further elucidate the conduction processes 
in this exciting novel type of molecular (super)conductor. 
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the "Stichting FOM" and is partially funded by the EC- 
RTN "QuEMolNa" (No. MRTN-CT-2003-504880), the 
EC-Network of Excellence "MAGMANet" (No. 515767- 
2) (Leiden), the DFG-Centre of Functional Nanostruc- 
tures (Karlsruhe), and the European Commission under 
the 6th Framework Programme through the Key Action: 
Strengthening the European Research Area, Research In- 
frastructures (No. RII3-CT-2004-505925). 
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